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ABSTRACT 



We report on the preliminary determination of the visual orbit of the double-lined 



<^ ' spectroscopic binary system BY Draconis with data obtained by the Palomar Testbed 

■ Interferometer in 1999. BY Dra is a nearly equal-mass double-lined binary system 



whose spectroscopic orbit is well known. We have estimated the visual orbit of 
. BY Dra from our interferometric visibility data fit both separately and in conjunction 

with archival radial velocity data. Our BY Dra orbit is in good agreement with the 
spectroscopic results. Due to the orbit's face-on orientation, the physical parameters 
_ implied by a combined fit to our visibility data and radial velocity data do not 

QQ ' yet result in precise component masses and a system distance, but with continued 

■ interferometric monitoring we hope to improve the mass estimates to better than 10% 

determinations . 

o 
o 
o 

1. Introduction 

Oh; 

Q ■ BY Draconis (HDE 234677, Gl 719) is a well-studied nearby (~ 15 pc), multiple stellar system 

containing at least three objects. The A and B components form a short-period (6 d) late-type 
. (K6 Ve ~ K7 Vvar) binary system, whose spectroscopic orbit is well known ( Bopp fc Evans 1973| , 



hereafter BE73; |Vogt fc Fekel 1979i hereafter VF79; and |Lucke fc Mayor 1980| , hereafter LM80). 



BY Dra is the prototype of a class of late-type flare stars characterized by photometric variability 
' due to star spots, rapid rotation, and Ca II H and K emission lines. Like the BY Dra system 

itself, a large fraction (> 85%, Bopp Sz Fekel 1977 , Bopp et al. 1980]) of BY Dra stars are known 



to be in short-period binary orbits. The rapid rotation of BY Dra A (period 3.83 d) that gives rise 
to the spotting and photometric variability is consistent with pseudosynchronus rotation with the 



A - B orbital motion (Hut 1981, Hall 1986), but pseudosynchronization is disputed by Glebocki & 
Stawikowski (1995, 1997) who assert asynchronous rotation and roughly 30° misalignment of the 
orbital/rotational angular momentum vectors. 

Despite the fact that the BY Dra A and B components are nearly equal mass, the system 
exhibits a significant brightness asymmetry in spectroscopic studies (VF79, LM80). VF79 
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attributes this to the hypothesis that A and B components of BY Dra are pre-main sequence 
objects (VF79, [Bopp et al. 19801 ), and are still in the contraction phase. VF79 argues for physical 
sizes of the A component in the range of 0.9 - 1.4 Rq, based primarily on rotation period and 
V sin i considerations. LM80 concur with the A component physical size argument from their v sin i 
measurements, estimating a 1.2 Rq size for a sinz ~ 0.5 (presuming rotation/orbit spin alignment 
with pseudosynchronization, and our orbital inclination from Table §). However, they continue by 
pointing out that if the A component macroturbulance were significantly larger than solar, then 
the V sin i measurements and the component diameters they are based on are biased high. If the 
pre-main sequence interpretation is correct, the BY Dra components are additionally interesting 
as an examples of the transition region between pre and zero-age main sequence states. 

BY Dra was detected as a hierachical triple system through common proper motion 
measurements of a BY Dra C component by Zuckerman et al. (1997); they find the C component 
separated by 17" from the A - B pair. Zuckerman's photometry on the C component is consistent 
with an M5 main-sequence interpretation {V - K ^ 6.2); assuming all three stars are coeval this 
clearly poses problems for the VF79 pre-main sequence hypothesis for the A and B components. 
At a projected physical separation of approximately 260 AU from the A ~ B binary, the putative 
low-mass C component would have negligible dynamical influence on the A - B binary motion. 
Further, the Hipparcos catalog ( ESA 1997|) implies the presence of at least one additional 



component as it lists BY Dra as having a circular photocentric orbital solution with 114 d period. 
This 114 d period is previously unreported in spectroscopic studies, and if correct it is difficult to 
understand why this motion was not previously detected. 

Herein we report on a preliminary determination of the BY Dra A - B system visual orbit 
from near- infrared, long-baseline interferometric measurements taken with the Palomar Testbed 
Interferometer (PTI). PTI is a 110-m H (1.6^m) and i^-band (2.2 fim) interferometer located 
at Palomar Observatory, and described in detail elsewhere ( polavita et al. 199S| ). PTI has a 
minimum fringe spacing of roughly 4 milliarcseconds (10~^ arcseconds, mas) in i^-band at the sky 
position of BY Dra, allowing resolution of the A - B binary system. 



2. Observations 

The interferometric observable used for these measurements is the fringe contrast or visibility 
(squared, V'^) of an observed brightness distribution on the sky. The reader is refered to one of 
the previous papers in this series (iBoden et al. 1999a| , |Boden et al. 1999b| , |Boden et al. 2000| ) for 



descriptions of the V"^ observables for binary stars. 

BY Dra was observed in conjunction with objects in our calibrator list (Table ||) by PTI in 
i^'-band (A ~ 2.2/um) on 21 nights between 23 June 1999 and 21 October 1999, covering roughly 20 
periods of the system. Additionally, BY Dra was observed by PTI in H-hand (A ~ 1.6/xm) on 13 
September, 2 October, and 20 October 1999. BY Dra, along with calibration objects, was observed 
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Object 


Spectral 


Star 


BY Dra 


Adopted Model 


Name 


Type 


Magnitude 


Separation 


Diameter (mas) 


HD 177196 


A7 V 


5.0 V/4.5 K 


6.6° 


0.70 ± 0.06 


HD 185395 


F4 V 


4.5 V/3.5 K 


9.9° 


0.84 ± 0.04 



Table 1: PTI BY Dra Calibration Objects Considered in our Analysis. The relevant parameters for 
our two calibration objects are summarized. The apparent diameter values are determined from 
effective temperature and bolometric flux estimates based on archival broad-band photometry, and 
visibility measurements with PTI. 



multiple times during each of these nights, and each observation, or scan, was approximately 130 
sec long. For each scan we computed a mean V'^ value from the scan data, and the error in the V'^ 
estimate from the rms internal scatter. BY Dra was always observed in combination with one or 
more calibration sources within ~ 10° on the sky. Table || lists the relevant physical parameters 
for the calibration objects. We have calibrated the data by methods discussed in Boden et 
al. (1998). Our 1999 observations of BY Dra result in 67 calibrated visibility measurements (50 
in K-hand, 17 in if-band). One notable aspect of our BY Dra observations is that its high 
declination (51°) relative to our Palomar site (33° latitude) puts it at the extreme Northern edge 
of the delay line range on our N-S baseline ( |Colavita et al. 19991 ), implying extremely limited 
u — V coverage on BY Dra. To our PTI visibilities we have added 44 double-lined radial velocity 
measurments: 14 from BE73, seven from VF79, and 23 CORAVEL measurements from LM80. 



Orbit Determination 



As in previous papers in this series ( Poden et al. 1999a , Boden et al. 1999b| , Boden et al. 2000) 



the estimation of the BY Dra visual orbit is made by fitting a Keplerian orbit model directly 
to the calibrated (narrow-band and synthetic wide-band) V"^ and RV data on BY Dra; because 
of the limited u — v coverage in our data derivation of intermediate separation vector models is 
impossible. The fit is non-linear in the Keplerian orbital elements, and is therefore performed 
by non-linear least-squares methods with a parallel exhaustive search strategy to determine the 
global minimum in the chi-squared manifold. The reader is refered to the previous papers for 
futher details on our orbit estimation procedures. 

Figure |l| depicts the relative visual orbit of the BY Dra system, with the primary component 
rendered at the origin, and the secondary component rendered at periastron. We have indicated 
the phase coverage of our V'^ data on the relative orbit with heavy lines; our data samples most 
phases of the orbit well, leading to a reliable orbit determination. The apparent inclination 
is very near the estimate given by VF79 based on the primary rotation period, assumed size, 
and assumption of parallel orbital/rotational angular momentum alignment. The orbit is seen 
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approximately 30° from a face-on perspective, which makes physical parameter determination 
difficult (§|). 

Figure ^ illustrates comparisons between our PTI V'^ and archival RV data and our orbit 
model. In Figure |^ six consecutive nights of X-band visibility data and visibility predictions 
from the best-fit model are shown, inset with fit residuals along the bottom. That there are only 
1-3 data points in each of the nights follows from the brief time each night that BY Dra is 
simultaneously within both delay and zenith angle limits. Figure |2|b gives the phased archival RV 
data and model predictions, inset with a histogram of RV fit residuals. The fit quality is consistent 
with previous PTI orbit analyses. 

Spectroscopic orbit parameters (from VF79 and LM80) and our visual and spectroscopic orbit 
parameters of the BY Dra system are summarized in Table ^. We give the results of separate 
fits to only our V'^ data (our 'V^-only Fit" solution), and a simultaneous fit to our V'^ data and 
the archival double-lined radial velocities - both with component diameters constrained as noted 
above. All uncertainties in parameters are quoted at the one sigma level. We see good statistical 
agreement between all the derived orbital parameters with the exception of the LM80 period 
estimate. 



4. Comparisons With Hipparcos Model 

The Hipparcos catalog lists a circular photocentric orbital solution for BY Dra with a 114 day 
period (ESA 1997), presumably in addition to the well-established 6 d period A - B motion. As 



noted above it is difficult to reconcile this hypothesis with the quality of the existing short-period 
spectroscopic orbit solutions from VF79 and LM80. However if the A - B system indeed did have 
a companion with this period, unlike BY Dra C it would lie within PTI's 1" primary beam, and if 
sufficiently luminous it would bias the visibility measurements used in our BY Dra A - B visual 
orbit model. We see no indications of this in our orbital solutions; the quality of our BY Dra 
visual orbit solution is consistent with our results on other systems. But it remains possible we 
have misinterpreted our V'^ data in the binary star model fit, and we are motivated to consider 
the 114-d periodicity hypothesis in the archival RV data. 

We note that the Hipparcos model is a photocentric orbit, and therefore calls for the A - B 
system to exhibit a reflex motion with radius > 0.05 AU at the putative distance of BY Dra. The 
114-d orbit hypothesis is at high inclination (113°), and therefore would produce a radial velocity 
semi-amplitude for the A - B system barycenter > 3.95 km s~^. This value is large compared to 
fit residuals observed by VF79 and LM80 in their spectroscopic orbital solutions (and by ourselves 
in the joint fit with our visibilities; Table §), suggesting the 114-d motion hypothesis is unlikely. 

To quantify this issue we have considered Lomb-Scargle periodogram analyses of the LM80 
BY Dra radial velocity data. We have chosen to use the LM80 data because it is the more precise 
sample, yielding an rms residual of roughly 0.5 km s~^ in our A - B orbit analysis. Figure S gives 
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Fig. 1. — Visual Orbit of BY Dra. The relative visual orbit model of BY Dra is shown, with the 
primary and secondary objects rendered at Tq (periastron). The heavy lines along the relative orbit 
indicate areas where we have orbital phase coverage in our PTI data (they are not separation vector 
estimates); our data samples most phases of the orbit well, leading to a reliable orbit determination. 
Component diameter values are estimated (see discussion in § ^), and are rendered to scale. 
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Fig. 2.— Joint yV^adial Velocity Fit of BY Dra. a: Six consecutive nights of V'^ data on 
BY Dra, and best-fit model predictions, b: Phased archival RV data (BE73, VF79, LM80) and RV 
predictions from our best-fit orbital model. Inset is a histogram of RV residuals to the fit model. 
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Orbital 
Parameter 


VF79 


LM80 


PTI 99 


y^-only Fit 


Full Fit 


Period (d) 


5.9750998 


5.975112 


5.975078 


5.975078 




± 8.7 X 10-5 


± 8.2 X 10-^ 




± 1.0 X 10-5 


To (MJD) 


41146.59 


43794.193 


51376.1407 


51376.1232 




± 0.07 


± 0.023 


± 0.0087 


± 0.0045 


e 


0.36 


0.3066 


0.321 


0.2978 




± 0.03 


± 0.0063 


± 0.015 


± 0.0032 


Ki (km s-i) 


28.2 ± 1.0 


28.55 ± 0.25 




28.31 ± 0.14 


K2 (km s-i) 


28.8 ± 1.8 


32.04 ± 0.35 




31.96 ± 0.15 


7 (km s""'^) 


-24.47 ± 0.65 


-25.35 ± 0.14 




-25.431 ± 0.078 


wi (deg) 


220 ± 5 


229.3 ± 1.3 


240.7 ± 8.2 


232.21 ± 0.91 


f^i (deg) 






106.7 ± 8.9 


115.1 ± 1.3 


i (deg) 






153.0 ± 3.6 


151.8 ± 3.5 


a (mas) 






4.376 ± 0.056 


4.399 ± 0.050 


AKcn 






0.69 ± 0.11 


0.571 ± 0.061 








1.39 ± 0.30 


0.92 ± 0.24 


Ay 




1.15 ± 0.1 






xVdof 






1.1 


2.4 (1.5 VV3-6 RV) 








0.033 


0.038 


\Rrv\ (kms-1) 




0.55 




2.3 (0.49 COR) 



Table 2: Orbital Parameters for BY Dra. Summarized here are the apparent orbital parameters 
for the BY Dra system as determined by VF79, LM80, and PTI. We give two separate fits to our 
data, with and without including archival double-lined radial velocities in the fit. Quantities given 
in italics are constrained to the Ustcd values in our model fits. We have quoted the longitude of 
the ascending node parameter ($7) as the angle between local East and the orbital line of nodes 
measured positive in the direction of local North. Due to the degeneracy in our V'^ observable 
there is a 180° ambiguity in O; by convention we quote it in the interval of [0:180). We quote mean 
absolute and RV residuals in the fits, \Ry2 \ and \Rrv\ respectively. 
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periodograms of the LM80 primary and secondary radial velocity data. First, Figure ^a gives a 
periodogram of the primary and secondary RV data, with probability of false alarm levels (P/a) as 
noted. As expected, both component lines exhibit a significant periodicity at the observed A - B 
orbit frequency of approximately (6 d)~^ (indicated by vertical line). In the same plot we sample 
the frequency of the 114-d orbit hypothesis, and no comparable periodicity is evident on the scale 
of the A - B motion. 

Presuming the 114-d motion hypothesis might be superimposed on the 6-d A - B orbit, in 
Figure ^b we give a periodogram of the LM80 primary and secondary RV fit residuals to the 
6-d hypothesis fit. In Figure ^3 we have adjusted to range of the periodogram to finely sample 
the frequency range around the (114 d)~^ hypothesis (indicated by vertical line). No significant 
periodicity is noted at this or any other frequency. The LM80 RV dataset spans roughly 400 days, 
so a 4 km s~^ amplitude periodicity in this dataset should have been evident in our analysis. 
Given these considerations, it seems unlikely that our V'^ measurements and A - B orbit model 
are affected by the presence of a third luminous body. 

5. Discussion 

The combination of the double-lined spectroscopic orbit and relative visual orbit allow us 
to estimate the BY Dra A - B component masses and system distance. However, because of 
the nearly face-on geometry of the A - B orbit, the accuracy of our inclination estimate, and 
consequently the component mass estimates and system distance estimate, is relatively poor. The 
low-inclination geometry is particularly difficult for astrometric studies because the astrometric 
observable {V'^ in this case) becomes highly insensitive to small changes in the inclination Euler 
angle. Table ^ lists the physical parameters we estimate from our Full-Fit orbit solution. The 
mass estimates are in a reasonable range for stars of this type but are formally crude - roughly 
24% 1-0" uncertainties, as is the system distance estimate - roughly 8% 1-a. In both cases the 
dominant error is the inclination; to improve these mass and distance estimates by a factor of two 
we must improve the inclination estimate by a factor of two. The orbital parallax with its large 
error is in 1-a agreement with the Hipparcos triginometric determination of 60.9 ± 0.75 mas, but 
we note that the Hipparcos solution was derived jointly with the 3 mas, 114-d orbital hypothesis 
that we believe to be suspect (§^. 

Table |3| also gives component absolute magnitudes and V - K color indices derived from 
archival broad-band photometry, our H and K component relative magnitudes, the LM80 V 
relative magnitude, and our system distance estimate. The component color indices are unaffected 
by errors in the system distance, and seem to be in good agreement with the color indices 
expected from stars in this mass range and the classical spectral typing of the BY Dra system. 
The absolute magnitudes are somewhat uncertain from the relatively poor distance estimate, but 
suggest a modest discrepancy with the mass-luminosity models of Baraffe et al (1998, BCAH98). 
As depicted in Figure ^, both components appear roughly 0.5 - 1 mag brighter at V and K than 
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Fig. 3. — Lomb Scargle Periodograms on LM80 BY Dra Radial Velocity Data. To assess the 
114-d orbit hypothesis for an additional companion to the BY Dra A - B system we have performed 
Lomb - Scargle periodogram analysis on double-lined BY Dra data from LM80. Top: standard 
periodograms on primary and secondary RV data from the LM80 dataset. Statistically significant 
response is seen at the well-known 6 d A - B period (vertical line). Bottom: Low-frequency 
periodogram of the residuals of the LM80 data to the best-fit A - B orbit model. No significant 
response is evident at the putative 114 d period (vertical line). 
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the BCAH98 models evaluated at our model masses would predict. The cause for this discrepancy 
could be abnormally large component sizes suggested by VF79, or the component masses are 
underestimated in our analysis. Clearly we need to improve our orbit model so as to further 
constrain the BY Dra component parameters. 

To assess the VF79 component size/pre-main sequence hypothesis, the most interesting 
measurement of the BY Dra A - B system would be unequivocal measurements of the component 
diameters. Unfortunately our data are as yet insufficient to determine these diameters 
independently. Canonical sizes of 1.0 Rq and 0.8 Rq (indicated by model effective temperatures 
and our IR flux ratios) yield angular diameters of approximately 0.6 and 0.5 mas for the A and 
B components respectively. At these sizes neither the H nor K-hand fringe spacings of PTI 
sufficiently resolve the components to independently determine component sizes. Consequently 
we have constrained our orbital solutions to these 0.6 and 0.5 mas model values. Our data does 
in fact prefer the slightly smaller primary component diameter to the larger 1.2 Rq size implied 
by VF79 and LM80 vs'mi and rotational period measurements. However either primary model 
diameter is possible with the expected systematic V'^ calibration errors. Additional data we will 
collect in the coming year (see below) may well place interesting upper limits on the component 
sizes, but unambiguous resolution of the BY Dra A and B components will have to wait for a 
longer baseline infrared interferometer; most likely the CHARA array currently under construction 
on Mt. Wilson Q 

The prospects for improving the BY Dra orbital model in general, and the orbital inclination 
estimate in particular are reasonable. To acheive 10% mass determinations on the BY Dra 
components we will need to reduce the inclination uncertainty to approximately 1.5°, which 
implies a rough quadrupling of the PTI visibility data on BY Dra, or an improvement of our V"^ 
measurement precision by a factor of two on a data set of similar size to that already collected. 
Both are plausible; PTI is expected to be in normal operation for at least another season, 
and hardware upgrades to the PTI fringe camera are planned and should provide significant 
improvements in sensitivity. 



The work described in this paper was performed at the Infrared Processing and Analysis 
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with the National Aeronautics and Space Administration. Interferometer data were obtained 
at Palomar Observatory using the NASA Palomar Testbed Interferometer, supported by NASA 
contracts to the Jet Propulsion Laboratory. Science operations with PTI are conducted through 
the efforts of the PTI Collaboration ( http: / /huey.jpl. nasa.gov / palomar / ptimembers.htm| ) , and we 
acknowledge the invaluable contributions of our PTI colleagues. We further acknowledge I.N. Reid 
(UPenn) for the compilation of mass-luminosity data on Henry & McCarthy and Andersen objects. 
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Fig. 4. — BY Dra Components in Mass/Luminosity Space. Here we give the positions of the A 
and B BY Dra components in observable Mass/My (upper pannel) and Mass/M^^ (lower pannel) 
spaces (an emulation of Figure 3 from BCAH98), superimposing low mass objects from Henry &: 
McCarthy (1993) and Andersen (1990) and BCAH98 solar-metalicity model tracks at ages of 10^, 
10^, and 10^ yrs. While our mass determinations are currently too crude to be definitive, both 
components appear over-luminous compared to the BCAH98 models, and it is difficult to reconcile 
the deviations with age effects. 
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Physical 


Primary 


Secondary 


Parameter 


Component 


Component 


a (10-2 AU) 


3.14 ± 0.36 


3.55 ± 0.41 


Mass (Mq) 


0.59 ± 0.14 


0.52 ± 0.13 


System Distance (pc) 


15.2 ± 1.2 


TTorb (mas) 


65.8 


± 5.4 


Model Diameter (mas) 


0.60 (± 0.06) 


0.50 {± 0.05) 


Mk (mag) 


4.47 ± 0.18 


5.04 ± 0.18 


Mh (mag) 


4.48 ± 0.21 


5.40 ± 0.29 


My (mag) 


7.48 ± 0.18 


8.63 ± 0.18 


V-K (mag) 


3.019 ± 0.039 


3.598 ± 0.061 



Table 3: Physical Parameters for BY Dra. Summarized here are the physical parameters for the 
BY Dra A - B system as derived primarily from the Full-Fit solution orbital parameters in Table 
0. Quantities listed in italics (i.e. the component diameters, see text discussion) are constrained to 
the listed values in our model fits. 



